subroutine disp_mol_info
use cell
use iounits
implicit none
   !----------------------------------------------------------------------------
   write( stdout, *  )
   write( stdout, 10 )
   write( stdout, 20 ) "Information on the selected Molecule: "//trim(molname)
   write( stdout, 30 ) "Molecule formula", trim(molform)
   write( stdout, *  )
   write( stdout, 40 ) "Total types of atoms", ntypeMol
   write( stdout, 40 ) "Total types of bonds", nbndtypMol
   write( stdout, 40 ) "Total types of angles", nangtypMol
   write( stdout, 40 ) "Total types of diherals", ndihtypMol
   write( stdout, *  )
   write( stdout, 40 ) "Total number of atoms", natomMol
   write( stdout, 40 ) "Total number of bonds", nbondsMol
   write( stdout, 40 ) "Total number of angles", nanglsMol
   write( stdout, 40 ) "Total number of diherals", ndihrsMol
   write( stdout, *  )
   write( stdout, 50 ) "Total weight of molecule", molmass
   write( stdout, 60 ) "Length of the molecule", minlen
   write( stdout, 70 ) "Total charge on molecule", sum(chargeMol)
   write( stdout, 80 ) "Experimental density", den0
   write( stdout, 10 )
10 format( 10x, 80("=") )
20 format( 10x, A       )
30 format( 10x, A30," :", 1x, A  )
40 format( 10x, A30," :", 1x, I5 )
50 format( 10x, A30," :", F10.4, " g/mol" )
60 format( 10x, A30," :", F10.4, " Angs"  )
70 format( 10x, A30," :", F10.4, " elec"  )
80 format( 10x, A30," :", F10.4, " g/cm3" )
return
end subroutine

subroutine disp_sys_info
use cell
use iounits
use wallinfo
implicit none
   !----------------------------------------------------------------------------
   write( stdout, *  )
   write( stdout, 10 )
   write( stdout, 20 ) "Information on the system: ", trim(molname), nmols
   write( stdout, 30 ) "Molecule formula", trim(molform)
   write( stdout, 40 ) "Total number of molecules", nmols
   write( stdout, *  )
   write( stdout, 40 ) "Total types of atoms", ntypeMol
   write( stdout, 40 ) "Total types of bonds", nbndtypMol
   write( stdout, 40 ) "Total types of angles", nangtypMol
   write( stdout, 40 ) "Total types of diherals", ndihtypMol
   write( stdout, *  )
   write( stdout, 40 ) "Total number of atoms", natom
   write( stdout, 40 ) "Total number of bonds", nbonds
   write( stdout, 40 ) "Total number of angles", nangls
   write( stdout, 40 ) "Total number of diherals", ndihrs
   write( stdout, *  )
   write( stdout, 50 ) "Total mass of the system", mass
   write( stdout, 60 ) "Total volume of the system", volume
   write( stdout, 65 ) "Box size in three dimensions", box
   write( stdout, 70 ) "Average density of the system", density
   if ( LAddWalls ) then
      write( stdout, 30 ) "Walls added on top and bottom", trim(WallName)
      write( stdout, 45 ) "X-Y dimension of the walls", nX, nY
      write( stdout, 46 ) "Layers on top and bottom", nZT, nZB
      write( stdout, 47 ) "Thickness of both walls", ThHig, ThLow
      write( stdout, 40 ) "Total atoms in the walls", sum(WallRemain)
   endif
   write( stdout, 10 )
10 format( 10x, 80("=") )
20 format( 10x, A, A, " x", I8   )
30 format( 10x, A30," :", 1x, A  )
40 format( 10x, A30," :", 1x, I8 )
45 format( 10x, A30," :", 1x, I3, " x", I3  )
46 format( 10x, A30," :", 1x, I3, 2x,   I3  )
47 format( 10x, A30," :", 1x, 2F10.4        )
50 format( 10x, A30," :", F15.4, " g/mol"   )
60 format( 10x, A30," :", F15.4, " Angs^3"  )
65 format( 10x, A30," :", 2(F8.2," x"), F8.2)
70 format( 10x, A30," :", F15.4, " g/cm^3"  )
return
end subroutine
